Fano-type interpretation of red shifts and red tails 
in hole array transmission spectra 
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We present a unifying point of view which allows to understand spectral features reported in 
recent experiments with two-dimensional arrays of subwavelength holes in metal films. We develop 
a Fano analysis of the related scattering problem by distinguishing two interfering contributions 
to the transmission process, namely a non-resonant contribution (direct scattering) and a resonant 
contribution (surface plasmon excitation). The introduction of a coupling strength between these 
two contributions naturally induces resonance shifts and asymmetry of profiles which satisfy simple 
scaling relations. We also report an experiment to confirm this analysis. 

PACS numbers: 78.20.Bh, 42.79.Dj, 73.20.Mf 
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The observation of extraordinary transmission of a 
metal film perforated with a two-dimensional periodic 
array of subwavelength holes 0,0 has been followed by 
a considerable amount of theoretical work jj, U, L* LA III • 
This work has not yet reached the maturity for a full 
understanding of the precise mechanisms involved in this 
transmission process. Important global characteristics of 
the spectra have not yet been addressed satisfactorily: 
we refer here to the fact that, in comparison with the 
naive "band structure" of surface plasmons (see equa- 
tion (J2J below) , resonances are red shifted and line shapes 
are asymmetric 0, 0, El- In this paper, we will address 
these issues by developing a unifying model along the 
lines of the original Fano analysis of the autoionization 
phenomenon in atomic physics [8| and by experimental 
confirmation of the essence of the model. We will show 
that Fano's treatment yields useful insight in the problem 
of transmission of light through metallic hole arrays by 
distinguishing two interfering contributions to the trans- 
mission process, namely a non-resonant contribution and 
a resonant one. The latter is associated with surface plas- 
mon (SP) excitations at the interfaces of the hole array. 
In order to simplify the physics, we allow for resonant SP 
excitation at one of the interfaces of the hole array only. 
This assumption is approximately valid for asymmetric 
arrays, i.e. arrays sandwiched between two different di- 
electric media (such as air and glass) in which a SP mode 
can be tuned into resonance on either the one or the other 
interface, but not on both at the same time. Fabry-Perot 
mediated coupling between two identical SP's living on 
separate interfaces can induce extra effects in the sym- 
metric case which we will not consider. 

The situation studied by Fano is one where the scat- 
tering from an input state can take place either directly 
towards a continuum of states (the scattered states) or 
via a quasibound state, that is through some type of reso- 
nant state, which is then coupled to the scattered states. 
This process therefore defines two distinct types of scat- 
tering channels: one open channel vb\ referring only to 



the continuum of states and one closed channel -02 with 
the resonant state, coupled to the open channel, as sym- 
bolized in figure A transition from the input state 
straight to the open channel ipi will be called 'direct' 
or 'non-resonant' as opposed to the other transition path 
that goes first through the quasibound state of the closed 
channel V>2 before being scattered and is therefore named 
'resonant'. The transition amplitudes associated with 
each path will interfere to define the total transition prob- 
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FIG. 1: a) Formal representation of the Fano model for cou- 
pled channels, b) Physical picture of the scattering process 
through the hole array directly (straight arrows) or via SP 
excitation. The curved arrows inside the hole represent the 
guided field. The period of the array is denoted do, its thick- 
ness h and the radius of each hole r. 

ability. This situation leads to typical asymmetric reso- 
nance line shapes, named (accordingly) as Fano profiles. 
We have to emphasize that Fano analysis only addresses 
the global characteristics of the line shape of the spec- 
trum through "structure parameters" but does not give 
theoretical expressions for these parameters, until a de- 
tailed model is defined. 

The process of transmission through metallic hole ar- 
rays is fundamentally a wave optics problem for which, in 
principle, a full electromagnetic approach could be given. 
However, this situation is, at the same time, a scattering 
problem which can be mapped onto a formal quantum 
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Hamiltonian problem, as soon as the various scattering 
channels arc identified. The scattering point of view re- 
lates the transition amplitudes that are derived from the 
Hamiltonian problem to the electromagnetic field scat- 
tering amplitudes. 

To calculate the total transition probability, one has 
to solve the coupled channel problem which is defined by 
the Hamiltonian 



H = H 1 +H 2 + V. 



(1) 



The eigenstates of H\ correspond to the continuum states 
of channel whereas H2 has one discrete state corre- 
sponding to the resonant state in channel ip2- The two 
channels are coupled via the coupling term V. Starting 
from an input state, the derivation of the eigenstates of 
the full Hamiltonian H will give us the transition am- 
plitudes from the input state to the new eigenstates of 
the coupled channel problem, taking into account the two 
possible paths. 

The first path is non-resonant and simply related to the 
direct scattering of the field through the subwavelength 
holes. This contribution will be referred to as Bethe's 
contribution, based on Bethe's solution to the problem of 
direct transmission of the electromagnetic field through 
an infinitely thin ideal metal with subwavelength holes 
II- We will simply adjust this solution to our situation 
of holes drilled on a real metallic film of finite thickness 
through a best-fit procedure, thus taking into account in 
a effective way the influence of these realistic conditions. 
This non-resonant path will correspond to the H\ Hamil- 
tonian formally describing a one-dimensional scattering 
problem, taken along the z direction, as shown in figure 
n Scattering states *f? E will be considered as eigenstates 
of H x , i.e. Hi^ p E = EV P E . They will be defined ac- 
cording to the direction p of the incident plane wave and 
are given asymptotically far away from the grating. The 
corresponding scattering amplitudes, related to Bethe's 
solution, can be gathered into a S- matrix Sb, which is 
symmetric (S B = Sb) due to transmission reciprocity 
property of the array. These states are non-resonant 
states (defined in the open channel) and belong as free 
fields to a two-folded continuum: one transmission con- 
tinuum and one reflection continuum. They are detected 
as a background in the reflection and transmission spec- 
tra. From the quantum point of view, these two continua, 
specific of an electromagnetic scattering process, are dis- 
tinguished by the quantum number p. 

The second path corresponds to a resonant contribu- 
tion, related to the existence of discrete SP excitations at 
the illuminated interface of the array. These excitations 
will be considered as discrete eigenstates ^ sp of the H2 
hamiltonian, i.e. H2^ S p = E sp ^ sp . It is well known that 
these excitations correspond to SP resonances defined on 
the periodic array by wavevector matching, making use 
of the array momentum wavevector [ic| . At normal inci- 



dence, such resonances are, to first order, given by com- 
plex resonance frequencies 
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where (n, m) are integers, ao is the period of the ar- 
ray and £1 and £2 are respectively the permittivities 
of the adjacent medium and the metal. The negative 
imaginary part of uj sp corresponds to the internal (non- 
radiative) damping of surface plasmon on a smooth in- 
terface. This damping represents the coupling between 
the surface plasmon and metal absorption losses. The 
real part of uj sp corresponds to the resonance frequencies. 
In practice, transmission spectra of hole arrays peak at 
frequencies that are typically 4% smaller than the value 
given by the real part of equation @ ■ This discrepancy, 
apparent in many experimental works 0. El 0- is well 
known though basically not understood so far. 

Due to the array periodicity, the incident field can be 
converted into a surface plasmon which will eventually be 
scattered (reflected or transmitted) by reciprocally cou- 
pling again with the array. This coupling between the 
discrete resonance and the scattering states is formally 
given by the V operator (see also figure QJd) which de- 
fines the radiative damping of plasmons at the surface 
of the array, considering the two scattering continua: a 
radiative transfer directly to free space (to the reflection 
continuum) and a radiative transfer through the holes 
(to the transmission continuum). One can give inter- 
pretations for these radiative dampings but the key re- 
sult of this paper relies only on the fundamental idea 
of the coupling with these two radiative channels. Any 
interpretation of this coupling, aiming at specifying the 
mechanisms of the radiative transfers, will not question 
this result; this is the inherent strength of a Fano-type 
analysis. 

Fano has solved the coupled channel problem giving 
the normalized eigenstate of the full Hamiltonian H 
[3. The coupling between the resonant state ^ sp (SP 
excitation) and the scattering state *& E , 
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is assumed hermitian. It introduces resonance shifts (A) 
and linewidths (r) whose expressions can be framed into 
a Kramers-Kronig relation 



dE'- 
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(4) 



The scattering dynamics of our system are fully con- 
tained in the S-matrix element corresponding to the 
global transmission amplitude 



t= ($E\S\i) 



(5) 



3 



with \i) denoting the input state, the left-incident field. 
This expression takes into account direct scattering from 
the input state to the continuum of scattering states ipi 
(Bethe contribution), specified through the S'-matrix Sb, 
as well as indirect scattering through the resonant state 
of the closed channel ip2- The expression for the S-matrix 
element (f5j) takes the shape of a sum of interfering terms 
from which the expression for the transmission coefficient 
(transition probability) immediately follows 



T=\t\ 2 = \t B 



a [E - (Esp + HA) + hS] 
[E - (E sp + HA)] 2 + [hT/2] 2 



(6) 



where E = huj. The 8 parameter defines the ratio be- 
tween the resonant transition amplitude and the direct 
(background) transition amplitude This expression 
for the global transmission coefficient reveals the super- 
position of the non-resonant 'Bethe' contribution, with 
the transmission coefficient |is| 2 , and of the resonant 
part, with a specific structure that turns out to be exactly 
a Fano-type profile when written in the natural variables 



E - (E sp + HA) _ 25 

hr/2 ' q ~ T' 



(7) 



The dimensionless q parameter determines the asymme- 
try of the profile and its sig n the corresponding blue 
(q > 0) or red (q < 0) tail [Tlf . 

Our analysis shows that the line shapes are fully deter- 
mined by 4 structure parameters: the resonance energy 
E sp , identified here as the real part of huj sp , the linewidth 
r, the resonance shift A and the asymmetry q. The 
linewidth T includes the plasmon damping terms, i.e. (i) 
radiative damping, defined by the coupling strength be- 
tween the surface plasmon discrete state and the scatter- 
ing states (reflection and transmission continua) , and (ii) 
non-radiative damping defined by the imaginary part of 
tvjj S p. The non-radiative damping contribution is gener- 
ally much weaker than the radiative one, so that the reso- 
nance shift A, given by the Kramers-Kronig relation 10}, 
is essentially due to the additional radiation damping. 
This shift can explain the discrepancy between expected 
resonances (see equation (JSJ) and actual peak positions. 
It defines therefore non-perturbatively the modified SP's 
dispersion relation when compared to the SP's disper- 
sion relation on a smooth interface, taking account of 
the "dressed" character of the SP by coupling to free-field 
modes. Note that, even if it can easily be parametrized 
in equations (|6I7|I . our one-resonance Fano description 
does not strictly allow for the inclusion of SP-SP scatter- 
ing since this would involve more than one discrete state. 
However, this contribution is of higher order in the limit 
taken here of small holes as compared to the wavelength. 

A full electromagnetic theory would deduce this new 
dispersion relation from the definition of the S'-matrix of 
the whole scattering process and from the extraction of 



its corresponding poles 0, 0] . The sign of the asym- 
metry q can be read off directly from the tail of an ex- 
perimental spectral profile of an isolated SP resonance. 
Finally, the line shapes of transmission profiles show a 
characteristic dip at the particular point e = —q corre- 
sponding to a destructive interference effect between the 
two channels. This destructive interference effect is con- 
comitant to the so called Wood anomaly [lij . 

The definition of the structure parameters is directly 
rooted in the introduction of the coupling between a non- 
resonant and a resonant channel, as essentially shown by 
the Kramers-Kronig relation This relation provides 
us with general relations between the damping and the 
resonance shift. Although we have no explicit expressions 
for the coupling strength Vg between surface plasmons 
and scattering states, we can argue generically how it de- 
pends on the (u),r,h) parameters. 

With a simple interpretative model, it is easy to see 
that the coupling strength is an increasing function of the 
frequency to. Towards the reflection continuum, the cou- 
pling is directly defined via the holes acting as Herz- or 
Rayleigh-type dipole emitters induced by the longitudi- 
nally polarized surface plasmon. Towards the transmis- 
sion continuum, the (evanescent) coupling through the 
holes is considered which is also an increasing function 
of the frequency. This monotonic behaviour of the cou- 
pling strength, that is of the linewidth T (see equation 
(|4"|)). implies, when inserted into this Kramers-Kronig re- 
lation that the resonance shift A is negative, i.e. a 
red shift, as noticed by several authors 0, 0,Hj| In fact, 
the absence of a detailed description for the linewidth- 
resonance shift ratio forces us to consider T as a free 
parameter, including both the effects of radiative damp- 
ing (surface plasmons coupled to scattering states) and 
non-radiative damping. 

This coupling is also strengthened when, at fixed fre- 
quency, the hole radius is increased. This lowers the 
cut-off frequency of the hole waveguide and actually in- 
creases the coupling of surface plasmons towards scatter- 
ing states of the reflection and transmission continua, as 
the radiating dipoles get stronger. The linewidth is there- 
fore predicted to be an increasing function of the hole 
radius (d r T > 0) and, at the same time, the red shift of 
the resonance should be more pronounced (d r \A\ > 0). 
These variations have been observed experimentally 0. 

If the thickness h of the array is increased, the cou- 
pling towards the transmission continuum decreases. In 
the hypothetic limit of infinite thickness, only radiation 
damping on the reflection interface remains. As h is in- 
creased, the progressive disappearance of the radiative 
channel through the holes, due to higher (non-radiative) 
losses inside the waveguide, induces a narrowing of the 
linewidth (dhT < 0) and consequently a reduced red shift 
(dh \A\ < 0) of resonance. Two regimes of variation for 
linewidths and resonance shifts can be understood from 
the introduction of a typical thickness above which the 



4 



radiative coupling is mainly towards the reflection con- 
tinuum. The variations dhT < and dh |A| < 0, along 
with the existence of a critical thickness, are again in 
agreement with recent measurements |l5). 

Figure [21 shows measurement transmission data from 
one of our experiments, for the (1, 0) SP resonance on the 
air-metal side of our array. This array is made of a 200 
nm thick gold (Au) film perforated with a square lattice 
of 70 nm radius holes spaced with a 700 nm lattice period. 
Our sample is asymmetric: the Au film is evaporated on 
a glass substrate with a 10 nm thick bonding layer of 
titanium (Ti). Strong absorption in this bonding layer 
prevents any surface plasmon from being excited on the 
metal-glass interface and, as such, guarantees that only 
one interface contributes, in agreement with the general 
frame of our model. This profile has a wavelength red 
tail, i.e. q is negative [l]]. We notice that the position of 
the measured peak (at A = 739 nm) is ~ 4% red shifted 
when compared to the position of the corresponding SP 
resonance A sp = 712 nm naively expected from equation 
(J2J with a period of the array clq = 700 nm. 




wavelength [nm] 



FIG. 2: Experimental transmission spectrum of the air-metal 
(1,0) SP resonance, obtained with an Au film of thickness 
200 nm with ao = 700 nm and r = 70 nm. The smooth 
curve shows the fitted Fano profile JSJ by adjusting freely 
the natural variables of Q. The vertical dotted line marks 
the position of the peak whereas the arrow shows the posi- 
tion of the resonance given by equation JS). The horizontal 
dotted line shows the Fano background level. Inset: Experi- 
mental transmission spectrum for the same array, but now in 
a film only 100 nm thick. The former dotted "cross" has been 
superimposed, emphasizing the background increase and the 
dispersive line shape. 

We attribute this shift to a negative A parameter, as 
to be expected from the above general arguments. It 
should be noticed that the various parameters of equa- 
tion © are intertwined in their effects on the line shape. 
Shifts of SP resonances are determined for one part by 
A and for the other part by the asymmetry itself. In 



natural units, the peak of the Fano profile is centered on 
e = l/q. A red tail asymmetry (q < 0) induces there- 
fore an extra red shift, though small (less than 1%) when 
compared to the red shift that is expected from A. On 
the other hand, the non-resonant transmission coefficient 
scales as |is| 2 ~ (r/A) 4 , according to Bethe's theory |IJ. 
This defines a non-resonant background colored towards 
the blue. The Fano peak will therefore be blue-shifted 
when superposed on this background. This blue shift 
can again be neglected when compared to the red shift 
A. 

In figure we have superposed a fitting of the experi- 
mental curve with a profile defined by equation JHJ 16] . 
The fitted parameters indicate a red shift of 4% deter- 
mined by A < and a relative linewidth of about 3%, 
that is of the order of the shift. The background is found 
to be only about 0.2%, which is consistent with Bethe's 
calculations for this system. Below 710 nm, the tail of 
the air- metal (1,1) resonance is emerging. 

As mentioned above, the value of the Fano asymmetry 
q characterizes the ratio between the resonant transition 
amplitude and the background transition amplitude. Our 
experimental curve confirms this scaling with an asym- 
metry of the order of 4 which can also be deduced from 
the fitting parameters. As a check of the Fano descrip- 
tion, we have reduced the Au film thickness from 200 nm 
to 100 nm with the same 10 nm thick Ti layer. With this 
reduction, the film thickness is only a few times the Au 
skin depth allowing another direct leakage channel. This 
should increase the non-resonant background contribu- 
tion and thus lower the value of q. We have confirmed 
this experimentally: see the inset of figure [21 where the 
background level is increased and the profile of the reso- 
nance gets close to a dispersive shape which corresponds 
to a lower q value. 

In conclusion, we have shown that a Fano analysis is a 
unifying tool to understand hole array spectra since sim- 
ple scaling laws can be inferred. Only one resonance has 
been addressed in this paper and one should now extend 
this work to the collection of SP excitations on the array, 
i.e. to multichannel theory. 

Note added in proof - During the process of publication 
of this paper, we learned about the theoretical work of 
Sarrazin el at 17] which shows important overlapping 
with our. Nevertheless, we put here the emphasis on the 
concept of coupled scattering channels in the context of 
experimental spectra, in order to understand correspond- 
ing lincshapes. 
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